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e$ • Abstract 

In order to answer this question, we combine ten independent astrophysical constraints 
in the space of the density parameters Q m of gravitating matter and f^A of vacuum energy. 
! We find that Q m = 0.31±0.07, tt A = 0.63±0.21, and thus fi m + fi A = 0.94±0.22. The to- 

tal x 2 is 4.1 for 8 degrees of freedom, testifying that the various systematic errors included 
are generous. We also determine VL m in the exactly flat case. Five supplementary flat-case 
' constraints can then be included in our fit, with the result Q m = 1 — f^A = 0.337 ± 0.031. 

§ ' It follows that the age of the Universe is t = 13.5 ± 1.3 (0.68/h) Gyr. 

O : 
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1 Introduction 

If the dynamical parameters describing the cosmic expansion were known to good pre- 
cision, we would know whether the Universe is open or closed, or whether its geometry 
is in fact exactly flat as inflationary theory wants it. To know the answer we need at 
least (i) the Hubble constant Hq, usually given in the form H = 100/i km s" 1 Mpc" 1 , 
(ii) the dimensionless density parameter Q m of gravitating matter, comprising baryons, 
neutrinos and some yet unknown kinds of dark matter, and (iii) the density parameter 
°f vacuum energy, related to the cosmological constant A by 

tt A = A/3H 2 . (1) 

A flat universe is defined by the condition 

n m + n A = i . (2) 

When Ho,Q m and Q A are known, the age of the Universe, to, can be obtained from 
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the Friedman-Lematre model as 

1 r 1 

t Q = — dx (1 - Q m - n A ) 
-Ho jo 

+n mX - 1 + n A x 2 r 1/2 - (3) 

In a previous publication (Roos & Harun-or-Rashid 1998) we tried to determine the 
preferred region in the (£l m , n A )-plane by combining three independent observational con- 
straints and a value for H . Over the years more observational constraints have become 
available, (some of them summarized in our unpublished preprint Roos & Harun-or-Rashid 
1999) so that we now can make use of fifteen independent constraints meeting our criteria. 
Since we combine the data in a least-squares fit, we can only make use of observations 
quoting a value and an error, but in addition many interesting limits also exist. 

In Section 2 we describe the fifteen observational constraints entering our least-squares 
fit. Ten constraints are valid in the space of Q m and Q A ; the remaining five are only 
valid on the flat line Eq. @ and will be included only when we fix the fit to that line. In 
Section 3 we describe the results of our fit in the (fl m , ^A)-pla,ne as well as along the flat 
line Eq. (0). We then also use Eq. fl3|) to determine t . 

2 Observational Constraints 

For the Hubble constant we use the value h = 0.68±0.05 from the analysis of Nevalainen & 
Roos (1998) in which the Cepheid period-luminosity relation is corrected for metallicity 
dependence. This agrees well with more recent precise determinations. For instance 
Mould & al (1999) find, when similarly corrected, h = 0.67 (our evaluation). 

2.1 Cosmic Microwave Background Radiation 

The observations of anisotropics in the CMBR are commonly presented as plots of the 
multipole moments Ct against the multipole £, or equivalently, against the FWHM value 
of the angular anisotropy signal. In general, the theoretical models for the power spectrum 
may depend on up to 9 parameters. Lineweaver (1998) and Tegmark (1998) have combined 
the data from MAP and PLANCK into a confidence region in the marginal subspace of 
the (Q m , f2 A )-plane. An independent analysis, combining the angular power spectrum of 
the BOOMERANG experiment (Melchiorri et al. 1999) with that of COBE, furnishes us 
a second constraint. 
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2.2 Gas fraction in X-ray clusters 

Matter in an idealized, spherically symmetric cluster is taken to be made up of a nearly 
hydrostatic inner body surrounded by an outer envelope, infalling with the cosmic mix of 
the components. The baryonic component of the mass in galaxy clusters is dominated by 
gas which can be observed by its X-ray emission. Thus by measuring the gas fraction, one 
expects to obtain fairly unbiased information on the ratio of Q m to the cosmic baryonic 
density parameter £l b . Using a very large sample of clusters, Evrard (1997) has obtained 
a 'realistic' value of 

^h-^^ (11.8 ±0.7) . (4) 

Taking Q b = 0.024 ± 0.006/i -2 from the low primordial deuterium abundance (Tytler, 
Fan & Buries 1996), one obtains 

tt m = 0.36 ±0.09 , (5) 

which we use as our constraint. 

This constraint are restricted to a flat Universe and we only use them together with 
the assumption of flat cosmology in section 3. 



2.3 Cluster mass function and the Lya forest 

In theories of structure formation based on gravitational instability and Gaussian initial 
fluctuations, massive galaxy clusters can form either by the collapse of large volumes in 
a low density universe, or by the collapse of smaller volumes in a high density universe. 
This is expressed by the cluster mass function which constrains a combination of Q m and 
the amplitude a of mass fluctuations (normalized inside some volume). The amplitude 
a is given by an integral over the mass power spectrum. From an analysis by Weinberg 
& al. (1998) combining the cluster mass function constraint with the linear mass power 
spectrum determined from Lya data, one obtains the relation 

fi m ±0.18Ov = 0.46 ±0.08 (6) 

which we use as one constraint. 
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2.4 X-ray cluster evolution 

Clusters of galaxies are the largest known gravitationally bound structures in the Universe. 
Since they are thought to be formed by contraction from density fluctuations in an initially 
fairly homogeneous Universe, their distribution in redshift and their density spectrum 
as seen in their X-ray emission gives precious information about their formation and 
evolution with time. Thus by combining the evolution in abundance of X-ray clusters 
with their luminosity-temperature correlation, one obtains a powerful test of the mean 
density of the Universe. 

The results of Bahcall, Fan and Cen (1997) can be summarized in the relation 

tt m = 0.195 ±0.11 + 0.071ft A (7) 

which we use as one constraint. 

The results of Eke & al. (1998) can be summarized in the relation 

tt m = 0.44 ± 0.20 - 0.077fi A (8) 

which we use as one constraint. 

Donahue & Voit (1999) constrain Vl m through a maximum likelihood analysis of te- 
peratures and redshifts of the high redshift clusters from the Extended Medium Redshift 
Survey, as well as from a low redshift sample (Markevitch 1998), finding 

tt m = 0.27 ±0.10 , (9) 

This constraint are restricted to a flat Universe and we only use them together with the 
assumption of flat cosmology in section 3. 

2.5 Gravitational lensing 

The number of multiply imaged QSOs found in lens surveys is sensitive to f2 A - Models of 
gravitational lensing must, however, explain not only the observed probability of lensing, 
but also the relative probability of showing a specific image separation. The image sepa- 
ration increases with increasing a*, the characteristic velocity dispersion. Thus the results 
can be expressed as likelihood contour plots in the two-dimensional parameter space of 
a* and Vt m = 1 - tt A (Chiba & Yoshii 1999). 

We integrate out a*, and we thus obtain the constraint 

tt A = 0.70 ±0.16 (10) 
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Im, Griffiths & Ratnatunga (1997) use seven field elliptical galaxies to determine 

tt A = 0.64 ±0.15 (11) 



in the flat model. 

We use these two constraints together with the assumption of a flat cosmology in 
section 3. 

2.6 Classical double radio sources 

There are two independent measures of the average size of a radio source, where size 
implies the separation of two hot spots: the average size of similar sources at the same 
redshift, and the product of the average rate of growth of the source and the total time 
for which the highly collimated outflows of that source are powered by the AGN. This 
outflow leads to the large scale radio emission. The two measures depend on the angular 
size distance to the source in different ways, so equating them allows a determination of 
the coordinate distance to the source which, in turn, can be used to determine pairs of 
Q m , Q a- values. 

By using 14 classical double radio galaxies, Daly, Guerra & Wan Lin (1998) determine 
an approximately elliptical 68% confidence region in the (fl m , ^A)-plane centered at (0.05, 
0.32). In our fit we use this constraint. 

2.7 Supernovae of type la 

Type la supernovae can be calibrated as standard candles, and have enormous luminosi- 
ties. These two features make them a near-ideal tool for studying the luminosity-redshift 
relationship at cosmological distances. 

The factor relating brightness to redshift is a function of Q m and Q\. The High-z 
Supernova Search Team Riess et al. (1998) have used 10 SNe la in the redshift range 0.16 
- 0.62 to place constraints on the dynamical parameters and to- 111 the (fl m , ^A)-plane 
their 68% confidence region is an ellipse centered at (0.20, 0.65) which we use as one 
constraint. 

The Supernova Cosmology Project (Perlmutter et al. 1998) has published an analysis 
based on 42 supernovae in the high-redshift range 0.18 - 0.83. In the (fl m , ^A)-pl& n e the 
68% confidence range is an ellipse centered at (0.75, 1.36) which we use as one constraint. 
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2.8 Power-spectrum of extragalactic objects 

Matter in every direction appear to be distributed in high-density peaks sepatated by 
voids. The average separation distance is ~ 130/i _1 Mpc, which translates into a peak in 
the power spectrum of mass fluctuations. This provides a co-moving scale for measuring 
cosmological curvature. Broadhurst & Jaffe (1999) used a set of Lyman galaxies at z ~ 3 
finding a constraint of the form 

tt m = 0.20 ±0.10 + 0.34fi A (12) 

Roukema & Mamon (1999) have carried out a similar analysis of quasars, finding 

tt m = 0.24 ±0.15 + (0.10 ±0.08)fi A (13) 

We use these two results as constraints. 



2.9 Galaxy peculiar velocities 

The large-scale peculiar velocities of galaxies correspond via gravity to mass density fluc- 
tuations about the mean, and depend also on the mean density itself. Two catalogs of 
galaxies have been analyzed for these velocities in order to provide information on Q m : 
the Mark III catalog (Willick & al. 1997) of about 3000 galaxies within a distance of 
~ 70k" 1 Mpc, and the SFI catalog (Borgani & al. 1999) of about 1300 spiral galaxies in 
a similar volume. Combining the results in these catalogs, Zehavi & Dekel (1999) quote 
the constraint 

tt m hl£n 2 ~ 0.58 ±0.12 , (14) 

in the case of flat cosmology, where the error corresponds to a 90% confidence range. 
Taking the index n of the mass-density fluctuation power spectrum to be n = 1.0 ±0.1 
(Bond & Jaffe 1998), one obtains the constraint 

tt m = 0.55 ±0.14 , (15) 

where the error corresponds to a 68% confidence range. This constraint we only use 
together with the assumption of a flat cosmology in section 3. 
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3 Results and Discussion 



3.1 Fits 

We perform a least-squares fit to the above ten constraints in the space of the two free 
parameters fl m , fl\- The Hubble constant is not treated as a free parameter, but is fixed 
to the value Nevalainen & Roos (1998). 

Paying rigorous attention to statistical detail, we use the standard minimization pro- 
gram MINUIT (James & Roos 1975). The best fit value is then found to be 

tt m = 0.31 ±0.07, A = 0.63 ± 0.21, 

X 2 = 4-1 • (16) 

In Fig.l we plot the shape of the la and 2a contours. 
The above values can be added to yield 

tt m + tt A = 0.94 ± 0.22 . (17) 

From this we conclude that (i) the data require a flat cosmology, (ii) the Einstein-de 
Sitter model is very convincingly ruled out, and (iii) also any low-density model with 
Q\ = is ruled out. 

If we assume exact flatness and refit the previous ten constraints as well as the five 
one-dimensional constraints (5), (6), (10), (11), (15), the result is 

tt m = 0.337 ±0.031, tt A = 0.663 ±0.031, 

X 2 = 7.3 . (18) 

The value of fl\ in the two-dimensional fit, Eq. (|i~6"D , is determined mainly by the 
constraint from the Supernovae Cosmology Project (Perlmutter et al. 1998). However, 
on the flat line several other constraints contribute much more strongly, so that the result 
in Eq. (|i~8f ) is very precise, even indepently of the supernovae constraint. 

Let us now substitute the above parameter values into Eq. (|). Adding a 7.4% H 
error quadratically to the density parameter errors in Eq. (|ljj), propagated through the 
integral (|3|), we find as a value for the age of the Universe 

to = 13.5 ± 1.3 (O.68//1) Gyr . (19) 

For an exactly flat Universe, only the error changes slightly to 1.1. 
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3.2 Systematic errors 

With results as precise as those for the flat model, a question arising is, what about 
neglected systematic errors? 

The constraints we use are indeed pulling the results in every direction. CMBR is 
orthogonal to the supernova constraints, the gravitational lensing constraint is exactly 
orthogonal to the gas fraction in X-ray clusters and the remaining constraints represent 
bands of several different directions. Thus we think that it is justified to consider the 
systematic errors as random and the total effect of possibly neglected systematic errors 
to be mutual cancellation. 

Moreover, the two-dimensional fit, Eq. (|16"D , x 2 is 4.1 for 8 degrees of freedom, and 
in the one-dimensional fit, Eq. (0), x 2 is 7.3 for 14 degrees of freedom, much too low 
for statistically distributed data. Thus we can conclude that the various errors quoted 
for our fifteen constraints are not statistical: they have been blown up unreasonably by 
the systematic errors added, and there is no motivation for blowing them up further by 
adding arbitrary systematic errors. 



3.3 Comparison with other data 

There are some categories of data which we have not used, but to which it is nevertheless 
interesting to compare our results. 

Totani, Yoshii & Sato (1997) have tested cosmological models for the evolution of 
galaxies and star creation against the evolution of galaxy luminosity densities. They have 
found Q\ > 0.53 at 95% confidence in a flat universe. 

Falco, Kochanek & Munoz (1998) have determined the redshift distribution of 124 
radio sources and used it to derive a limit on Q m from the statistics of six gravitational 
lenses. In a flat universe their best fit yields Q m > 0.26 at 95.5% confidence. 

Since we also determine a value for the age of the Universe, it is of interest to look at 
other recent to determinations. By using several techniques Chaboyer (1998) determine 
the age of the oldest globular clusters (GC) to t GC = 11.5 ± 1.3 Gyr. Another estimate 
of the age of the oldest globular clusters is due to Jimenez (1998). He quotes the 99% 
confidence range tec — 13.25 ± 2.75 Gyr, which translates into the 68% confidence range 
tec — 13.3 ±1.1 Gyr. Taking a mid-value and adding 0.9 Gyr systematic error to account 
for the spread of values, we have 
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t GC = 12.4 ±1.3 ±0.9 Gyr . (20) 

In order to obtain to from the GC value, however, one must add the time it took for 
the metal-poor stars to form. Estimates are poor due to the lack of a good theory, so 
Chaboyer (1998) recommends adding 0.1 to 2 Gyr. We presume that this is a 90% CL 
estimate, +1 ± 1 Gyr. Thus one arrives at the 68% confidence range 

to = 13.4 ± 1.3 ± 0.9 ± 0.6 Gyr = 13.4 ± 1.7 Gyr . (21) 

Jimenez (1998) also finds an age of to = 13 ±2 Gyr for the galaxy 53W069 at z = 1.43. 
This value as well as Eq. (21) are in excellent agreement with our value in Eq. (|19[) 
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Figure 1: The la and 2a statistical confidence regions in the (fl m , f^-plane are shown. 
The '+' marks the best fit: (Q m ,QA) = (0.31,0.63). The diagonal line corresponds to a 
flat cosmology. 
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